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Abstract 

We study whether a violation of the null energy condition necessarily implies the 
presence of instabilities. We prove that this is the case in a large class of situations, 
including isotropic solids and fluids relevant for cosmology. On the other hand we 
present several counter-examples of consistent effective field theories possessing a stable 
background where the null energy condition is violated. Two necessary features of 
these counter-examples are the lack of isotropy of the background and the presence 
of superluminal modes. We argue that many of the properties of massive gravity can 
be understood by associating it to a solid at the edge of violating the null energy 
condition. We briefly analyze the difficulties of mimicking H > 0 in scalar tensor 
theories of gravity. 


1 Introduction 

Energy conditions play an important role in general relativity and are a basis of singularity 
theorems [1] and entropy bounds [2], These conditions reflect the microscopic properties of 
the medium sourcing the energy momentum tensor. In some cases they may be violated 
without making the system physically unacceptable. This is the case of the strong energy 
condition, whose violation in Nature is indeed suggested by the success of the theory of 
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inflation and also by the recent observation of cosmic acceleration. It is however widely 
believed that any well behaved physical systems should respect the weakest of these condi¬ 
tions stating that T^n^n 1 ' > 0 for any null vector n. This condition is often referred to as 
the null energy condition (NEC). For an isotropic fluid with positive energy density, NEC 
amounts to w = p/p > —1. In some simple cases a connection between violation of NEC 
and the presence of instabilities in the system sourcing T jW was already established (see for 
instance [3]). It would however be desirable to investigate under what condition this con¬ 
nection becomes a theorem. This question is obviously theoretically interesting 1 , but it may 
also turn out to be relevant to present day cosmology. Indeed, recent supernova data mildly 
favor a NEC-violating equation of state [5], w < — 1. If this result were to be confirmed 
by more precise experiments in the future, it would become essential to understand under 
what conditions NEC can be violated in the absence of pathologies. A no-go theorem valid 
in ordinary 4D theories may point towards more exotic scenarii, like for instance the Dvali- 
Gabadadze-Porrati model [6] where an effective w slightly less than —1 for dark energy can 
be achieved [7, 8], without manifest pathologies [9]. Finally, as we will show, the study of 
NEC violating field theories opens up a new perspective on theories of massive gravity (see 
Refs. [10, 11, 12, 13]). 

The purpose of the current paper is to study in detail the relation between instabilities 
and NEC violation. For most of the paper we will be dealing with a system of derivatively 
coupled scalar fields with coordinate dependent condensates. Although this system may 
seem somewhat artificial, its simplicity allows us to concentrate on the physically relevant 
features avoiding technicalities which arise when one includes higher spin fields (e.g. vectors). 
Moreover, as we will explain in section 5 our results are in fact significantly more general. In 
particular our arguments apply to arbitrary relativistic solids and fluids. The main reason 
for this generality is that NEC violation is related to instabilities in a particular sector of 
the theory. This sector consists of just the Goldstone degrees of freedom (sound waves) 
associated to the spontaneous breakdown of local space-time translations. The properties 
of the Goldstone modes are dictated by symmetry, and are therefore largely independent of 
the specific fields that gave rise to the background. 

The paper is organized as follows. In section 2 we describe the basic setup. In section 3 
we study the general conditions for stability. We concentrate on some unusual features of 
rotationally non-invariant backgrounds, which make it rather difficult to establish a rela¬ 
tion between NEC violation and instabilities in general. In particular, we explain that the 
positive definiteness of the Hamiltonian is not necessary for stability. Another peculiarity 
is related to the fact that in an effective field theory there is no a priori reason to exclude 
the presence of superluminal excitations around a background which spontaneously violates 
Lorentz invariance. It turns out that in the presence of superluminal modes stability is not 

1 For instance, NEC prevents the existence of Lorentzian wormholes [4] and of flat bouncing cosmologies. 
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a frame independent notion. More precisely a perfectly stable system may appear to have 
ghost or gradient instabilities in some Lorentzian reference frames. 

In section 4 we establish a relation between NEC violation and instabilities. We start by 
analyzing two-dimensional systems. We show that if Goldstone modes have non-degenerate 
dispersion relations 2 and if the energy-momentum tensor can be diagonalized by a Lorentz 
boost, then NEC violation implies the presence of an instability. On the other hand, we 
construct a simple 2-dimensional model which is stable (the Hamiltonian is even positive 
definite), violates NEC but possesses a non-diagonalizable energy momentum tensor. We 
then present generalizations of these results to higher-dimensional space-times. We prove 
that a violation of NEC leads to instabilities in a large class of situations. These include 
systems that are either rotationally invariant or possess only subluminal excitations. Never¬ 
theless we were able to construct an example of a stable low energy effective field theory that 
violates NEC. This system has a positive definite Hamiltonian. Unlike in the 2-dimensional 
example, the energy-momentum tensor of this system can be diagonalized by a Lorentz boost. 
A peculiar necessary feature of stable NEC violating systems is the presence of superluminal 
modes in the spectrum of perturbations. 

In section 5 we provide a more general perspective on our results. We relate scalar field 
theories with coordinate dependent condensates to relativistic solids and fluids, and elaborate 
on this correspondence, rederiving in field theory language some of the known results of 
hydrodynamics (such as Kelvin theorem). Although most of the results of this section are 
well known, we find it useful to present this discussion as it provides an intuitive physical 
picture of the systems considered in the previous sections and in works on massive gravity. 
A pragmatic reader can skip this section and proceed directly to section 6. In the latter we 
study NEC violation in rotationally invariant systems with degenerate dispersion relations. 
We argue that, unlike in massive gravity [13], these degeneracies cannot be protected by 
symmetries. 

In section 7 we speculate on the implication of our results were future more precise 
cosmological obervations to imply w < — 1. We present simple arguments explaining why it 
is difficult to mimic such an expansion in scalar-tensor theories of gravity without introducing 
a NEC violating energy-momentum tensor. This allows us to propose a relatively simple 
candidate model which can do the job without obviously contradicting the precise tests of 
general relativity. We summarize our results in section 8. 

It is worth noting that a recent paper [14] claims to have a proof that NEC violation im¬ 
plies instabilities in a broad class of systems including those considered in our work. However, 
we believe that proof is not correct as we have been able to provide explicit counterexamples 
that contradict it (see section 4). 

2 Degenerate dispersion relations, where some modes have either vanishing gradient term or vanishing 
time derivative term (or both), always emerge in the consistent theories of massive gravity [10, 11, 12, 13]. 
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2 Setup 

Let us consider a system of N scalar fields (j) 1 , i — 1 ,,N interacting with gravity 

S = Sgr + f d 4 x y/—g F^ 1 , B IJ ,...) , (1) 

where S gr is the gravitational part of the action, 

B IJ = g*” % (2) 


and the dots stand for higher derivative terms. We work in the context of effective field 
theory and assume that the scalar action has a cut-off A smaller than the Planck scale Mpi. 
More specifically we assume the function F in eq. (1) is of the form 




ij 


= A 4 F 


/ecj) 1 B IJ d 2 (j) 1 d 2 (p J 


V A ’ A 4 


A 6 


( 3 ) 


where e < 1 is a small parameter responsible for the “lightness” of the fields (j ) 1 , while all 
the other dimensionless parameters in the function F are of order one. Indeed for the very 
existence of a non trivial low energy effective field theory the mass of the scalars ~ eA should 
be much less than the cut-off A. Notice also that in the limit e —» 0 we would be dealing 
with a set of naturally massless Goldstone bosons. We want now to consider some classical 
solution , 0^ and study perturbations around it. The validity of our effective field theory, 
or the control over quantum corrections, forces us to consider only solutions that are smooth 
on space-time scales of order A -1 


C00 7 ] > 



( 4 ) 


In other words we assume that the fields do not vary appreciably over lengths of order 1/A or 
more formally ]<9 /t <A . Notice that this constraint does not necessarily imply 3 d^cj) 1 <C A 2 . 
Let us consider the relevant length scales. Since the energy momentum generated by the 
scalars is of order A 4 , space-time will be curved over a length tq rsj Mpi/K 2 . Under the 
assumptions A <Mp; and e<l there is a range of lengths 


1 

- -C r -C r G , 


1 

eA 


( 5 ) 


where we can use effective theory, spacetime is approximately flat and the scalars can be 
treated as massless. In physically interesting situations where the scalars describe some sort 

3 In the Goldstone limit, e = 0, the terms of the type [ B IJ ]" represent the leading interaction among n- 
quanta in the derivative expansion. In this respect \B IJ /A 4 \ is controlling the size of classical non-linearities 
rather than the validity of the quantum expansion, which is instead more directly tight to the condition 
| <90 <A (see a related discussion in Ref. [15]). 
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of cosmic fluid or a generalized quintessence, the above range would encompass all relevant 
sub-horizon scales. We want to study perturbations around the solution in the above range 
of scales (cf., e.g., Ref. [16]), that is for space momenta k 

-rj-t eA < k < A . (6) 

Mpi 

As we focus on a patch of space-time where curvature and scalar masses are negligible, by a 
suitable change of coordinates we can approximate by the Minkowski metric r] fll/ in the 
whole patch and the background scalar fields (ft^ by linear functions of the coordinates 

tii = 0o + K x>l • ( 7 ) 

Since we can neglect both masses and higher derivative terms we can define our function F 
to be just F = A 4 F(B IJ /A 4 ). Finally, at momenta larger than the curvature, k S> A 2 /Adpi, 
we can neglect the mixing between scalar fields and metric perturbations. Indeed this mixing 
is a manifestation of the gravitational Higgs mechanism: the mixing scalar modes are the 
Goldstone bosons of spontaneously broken translations, and the graviton has effectively a 
mass of the order of the curvature A 2 /Mpp By a suitable (unitary) choice of coordinates 
the Goldstones can be eliminated, in which case extra physical polarizations appear in the 
graviton held. However when working at k A 2 /A / Ipi, the graviton mass can be neglected 
and, by an analogue of the equivalence theorem, the extra graviton polarizations can be 
accurately described by the Goldstone bosons. In the rest of the paper we will work in A = 1 
units. 

By defining the perturbed scalar held as (ft 1 = (ft T d + n 1 we can write the quadratic 
Lagrangian for scalar perturbations in Fourier space as 


L = Lu (p) n 1 ( —p)'n" J (p) . 

Here p is the four-momentum and the matrix Ljj(p) is given by 

Lij = 2F IKiJL A^p^ A^p u + F LJ p 2 , 

where 


F u = 


d F 


d 2 F 


dB u ’ rj KL dB IJ dB KL ' 

It is straightforward to check that the background energy-momentum tensor T jW is given by 


( 8 ) 


(9) 


( 10 ) 


Tnu = ZFuA 1 Al - Fg 


( 11 ) 


Our goal is now to study the stability of eq. (8) when the background violates NEC, that is 
when there exists some light-like vector such that 


T^/rFrF < 0 . 


( 12 ) 
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3 General criteria for stability 

3.1 One flavor in two dimensions 

Before addressing the stability of Lagrangian (8) in full generality, it is instructive to discuss 
a simple special case: a single scalar field in (l+l)-dimensional space-time with action (fr = 
d t n, 7r' = d x ir) 



The dispersion relation is 

X(u) = Kv 2 + Mao — G = 0 , (14) 

where 

v = u>/\k\ (15) 

is the ratio of the frequency and the absolute value of the spatial momentum k, and a = ±1 
depending on the direction of the spatial momentum. 

In general, the system (13) can exhibit two different types of instability. First, there can 
be a gradient instability. This is present whenever equation (14) has complex roots. This 
is similar to the usual tachyonic instability which arises in a Lorentz invariant theory when 
some field has a negative mass squared, m 2 <0. It is however different since in the latter 
case the time-scale r with which the instability grows is set by the inverse mass, r ~ \rn\~ 1 . 
This time scale may be very long compared to the other characteristic time scales if m 2 
is small enough. Equivalently, conventional tachyonic instabilities are present only for long 
wavelength modes. On the other hand if eq. (14) has a complex root this implies the presence 
of instabilities at all wavelengths, so that their rate of growth is not bounded. Instabilities 
of this type are significantly more dangerous than the usual tachyonic ones. Nonetheless in 
the rest of the paper we will improperly refer to gradient instabilities as tachyons. 

Another type of instability is the one associated with ghosts. In a Lorentz invariant 
theory a ghost is a field entering the action with a wrong sign kinetic term. At the classical 
level such a field has a negative definite free Hamiltonian. Therefore, neglecting non linear 
terms or interactions with other fields, no runaway behaviour occurs: there just happens to 
exist one decoupled sector using the opposite sign convention for time and energy. Runaway 
behaviour shows up when the ghost is set to interact with ordinary fields having positive 
definite free Hamiltonian. Notice that, being due to non-linearity, the time scale for growth 
becomes infinite in the limit of vanishing initial perturbation, i.e. in the classical vacuum. 
The situation is generalized, but gets worst, at the quantum level, just because the held 
fluctuates away from the origin already in the vacuum. The ghost quanta have negative 
energy 4 , so that, once interactions are taken into account, the vacuum itself becomes unstable 
4 We quantize on a physical Hilbert space with positive norm, so that all probabilities are non-negative. 
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with respect to spontaneous creation of positive and negative energy quanta with vanishing 
total energy. Again, ghost instabilities are entirely due to the non-trivial couplings between 
negative and positive energy fields, and strictly speaking require to go beyond the quadratic 
approximation. We also recall that in a Lorentz invariant theory, the rate of development 
of such an instability is uncalculable (formally infinite), even for arbitrary small couplings, 
due to the non-compact final state phase space volume. A way to solve this problem is 
to introduce a cutoff which explicitly breaks Lorentz invariance. In such a way the rate 
of growth of the instability can be at least estimated and phenomenological consistency 
satisfied if the coupling between the ghosts and the other fields is small enough (see for 
instance refs. [17, 18]). In what follows we will not be concerned with the latter possibility 
but instead require the stronger, more robust, condition that ghost instabilities be absent. 
As made clear from the discussion, whether or not negative energy states in the system 
(13) imply instability depends on the spectrum of the fields coupled to them. One coupling 
that surely exists is the one to the graviton. Moreover in the range we are focusing on, the 
graviton satisfies to a good approximation a Lorentz invariant dispersion relation u 2 = k 2 . 
We will therefore require the system to be kinematically stable with respect to the associated 
production of tt quanta and of some other massless particle (in particular the graviton) 5 . 

Let us now study what the stability conditions are. Absence of tachyons requires that 
both roots of the equation (14) be real, i.e. 

M 2 + 4 KG > 0 . (16) 

The issue of ghosts can be equivalently studied at the classical or quantum level. We find 
it convenient to work right away at the quantum level and use the perhaps more intuitive 
particle physics viewpoint. The general solution of the held equations following from the 
action (13) has the form 

7 r(t,x) = I dk (ai(Jfe) e ik{uit ~ x) + a 2 (k) e ik{u2t ~ x) ) , (17) 

where z/ 12 are the two solutions of eq. (14) with a — 1. Without loss of generality, we assume 
u\ > iJ- 2 - Reality of the held tt implies 

ai,2(ky = ai^(-k) . (18) 

By eq. (13) the momentum II canonically conjugated to tt is given by 

n = Kir + \Mtt' . (19) 

5 Note, that in two dimensions the graviton does not have propagating degrees of freedom. This is not 
essential for our purposes as we will eventually be interested in four dimensional theories. In the 2d example 
we are therefore considering the coupling to some other massless field. 
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The canonical commutation relations in momentum space become 


a a (k),a b (k')\ 


25 ab a(K)(-l) a 
7T ky/M 2 + 4 KG 


S(k + k') , 


( 20 ) 


where a(K ) stands for the sign of K. It is convenient to consider separately the two cases 
K > 0 and K < 0. 


• K > 0 . 

In this case we define creation and annihilation operators a'(k), a(k ) as 

a\k) = a\(k)6(k) + oi 2 (k) 0 (—k) (21) 

a{k) = a>i(—k)9(k) + a 2 (— k)6(—k) (22) 

Then, by defining the Fock vacuum as the state |0) satisfying 

a(/c)|0) = 0 Vfc, (23) 

the Fock space generated by the action of the a'(/c)’s on |0) has a positive definite inner 
product. Consequently, the energy spectrum of one particle states is 

E{k) = V\k6{k) + u 2 k9(—k) . (24) 

Negative energy states are absent provided U\ is positive and z/ 2 is negative. For K > 0 this 
condition is true if and only if G > 0. Notice that positivity of K and G coincides with 
positivity of the classical Hamiltonian H 

H = I dx^ (K-\ n - ±Mtt') 2 + Gtt' 2 ) > 0 . (25) 

However, K,G > 0 is not necessary for the stability of the system (cf. Ref. [15]). Indeed, 
consider the situation where V\ and are both positive, in which eq. (24) is not positive 
definite. The ghost instability is still absent, provided the two-momentum of the multiparticle 
states of the 7T-system (13) never belongs to the interior of the lower light cone, that is if 

E>-\k\. (26) 

In this case, in spite of the presence of negative energy states, vacuum production of negative 
energy 7r-particles and positive energy gravitons is not kinematically allowed. The basic 
reason for stability in this case is that in addition to energy there exists an extra conserved 
charge, the total momentum. Condition (26) holds even when both iq and v 2 are positive 
provided 


v-i < 1 , 


( 27 ) 



or if both v\ and V 2 are negative provided 

z/i > -1 . (28) 

In both cases the Hamiltonian is unbounded both from above and from below. 

To better understand the physics of such systems let us consider a subluminal scalar held 
in the reference frame where its action is rotationally invariant 6 , 

S — Jdtdx ^ ( 7 I 2 — v 2 7i' 2 ) , (29) 

with 0 < v 2 < 1. The Hamiltonian associated to this action is positive definite and there is 
no doubt that this system is perfectly stable. Let us, however, consider the same system as 
seen in another Lorentzian reference frame, i.e. let us perform a change of time and space 
variables in the action (29) according to the rule 7 

0M£ sx:)- 

where ch = cosh A and sh = sinh A are the boost parameters. After this change of coordinates 
the action (29) takes the general form (13) with the coefficients K, M, and G equal to 

K = ch 2 — sh 2 v 2 

M — | ch sh (v 2 — 1) (31) 

G = ch 2 v 2 — sh 2 

One observes that for v 2 < 1 the coefficient K remains positive in all reference frames, 
while G becomes negative for large enough boosts, sh 2 /ch 2 > v 2 . So, in the boosted frame, 
negative energy states appear and the Hamiltonian is unbounded both from above and from 
below. Dispersion relations for one-particle states in the original and in the boosted reference 
frames are shown in Fig. 1. It is clear from this figure that the conditions (27), (28) hold in 
both frames as it is natural to expect for a stable system. 

Actually, the presence of negative energy states is easy to understand in the following 
way. Let us consider a massive particle with relativistic dispersion relation 

E 2 = k 2 + m 2 , (32) 

moving with velocity V in the original reference frame. At V > v the emission of Cerenkov 
radiation is kinematically allowed and the massive particle can emit one quantum of the n- 
ficld, thus losing part of its energy. On the other hand, from the point of view of an inertial 
6 In two dimensions, rotations reduce to the Z 2 reflection x —> —x. 

7 Recall that all our considerations are in the context specified in section 2: starting from a Lorentz 
invariant action, we study the dynamics of fluctuations around a generic Lorentz-breaking field configuration. 
Lorentz transformations are still a symmetry of the whole system. 
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Figure 1: The dispersion relations (bold lines) of a subluminal scalar in the original frame 
in which eq. (29) holds (left), and in a highly boosted frame (right). 

observer moving at velocity V, the massive particle is initially at rest and then starts moving 
by emitting a negative energy quantum. 


• K < 0 . 

Let us now come back to the study of the general action (13) and consider the case K < 0. 
It is straightforward to check that the energy spectrum of one-particle states is now given 
by 

E(k) = u\k 9(—k) + u 2 k 9(k), (33) 

so that a state with two particles of equal and opposite spatial momentum has negative total 
energy (recall that zq > zq). At first sight, this seems to necessarily imply that theories with 
K < 0 suffer from a fatal ghost instability. 

However, this naive conclusion is not quite correct. Indeed, let us again consider action 
in eq. (29) but for the case of a superluminal 7r (v 2 > 1). There is no doubt that this action 
describes a perfectly stable system in this case as well. Let us now look at this system in 
the boosted frame (30). eqs. (31) now imply that the kinetic coefficient K becomes negative 
for large enough boost parameters, ch 2 /sh 2 > v 2 (see Fig. 2). 

Naively, this result leads to a paradox where the vacuum of a stable system (29) seems 
unstable to a fast moving observer. Let us assume for the sake of the argument that there 
exists a vertex with two gravitons and one n. Then, compatibly with conservation of 2- 
momentum, the moving observer expects to see events of spontaneous creation of one negative 
energy 7r-particle and a pair of positive energy gravitons from the vacuum, see Fig. 3. He 
should then conclude that the vacuum is unstable. 

The resolution of the paradox is just that the observer at rest and the moving one 
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Figure 2: The dispersion relations (bold lines) of a superluminal scalar in the original frame 
(left) and in a highly boosted frame (right). 

have performed mutually inequivalent quantizations of the system. This is understood by 
considering the causal structure of the system and the initial conditions associated to the two 
different parametrizations of time. First notice that, for a large enough boost, part of the 
past causal cone is mapped to the future of the new space-time coordinates. Consequently, 
an incoming 7r-particle in the original (unboosted) reference frame, is mapped to an outgoing 
7r-ghost in the boosted frame. Therefore the decay process n —» 2g becomes the vacuum 
decay process vacuum ^ it + 2g in the boosted coordinates. Now the point is that the 
two observers have defined causality by two inequivalent Cauchy surfaces, corresponding to 




Figure 3: The creation of a negative energy particle and two gravitons out of vacuum in the 
boosted frame (right). The bold arrows are the particles momenta, which sum up to zero. 
The same process appears in the original frame as the decay of a superluminal particle into 
two gravitons (left). 
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Figure 4: The bold lines depict the past causal cones in the original reference frame (left) and 
in a highly boosted one (right). Left: The two distinct vacuum choices require no incoming 
wave in the light shaded region (t < 0) and in the dark shaded region ( t' < 0), respectively. 
As clear from the picture these choices are inequivalent: an incoming particle moving along 
the left branch of the past causal cone is discarded in the former choice while it is accepted 
(as outgoing) in the latter. 

surfaces of constant time in the two different coordinate systems. This is made evident in 
Fig. 4 where the two distinct vacuum choices correspond to the request of no-incoming wave 
in two distinct shaded regions. The two choices are inequivalent because there are incoming 
(outgoing) portions of particle trajectories that are allowed in one choice and not in the other. 
Another way to say it is that the time vector, which defines the hamiltonian flow (and its 
quantization), is inequivalently chosen in the two parametrizations. This inequivalence is 
reflected by the fact that the direction of time flow of one of the two particle trajectories 
that cross each spacetime point is mutually reversed in the two cases. As shown in Figure 
4 this is precisely like comparing two parametrizations of spacetime where the role of time 
and space is interchanged 8 . 

In the end, as it should, the behaviour of the system does not depend on the observer 
but on the choice of quantization procedure, i.e. on the vacuum state. Stability favors the 
obvious vacuum, the one associated to the time slicing operated by the observer at rest. 
Provided this is the vacuum, the moving observer will perhaps find it a little puzzling that 
he does not observe copious production of 7T+2 g from nothing. However, with some thinking, 
he will realize that because of his unhappy choice of time coordinate, the correct causal past 
overlaps with a space time region with infinitely growing time coordinate. 

8 Notice, by the way, that if we perform the exchange of time and space for an ordinary Lorentz invariant 
massive field in two dimensions, we find that it becomes a tachyon. 
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It is worth noting that there were suggestions [19] that boundary conditions at future 
infinity may help to solve the cosmological constant problem. It would be interesting to 
understand whether the situation just discussed may be relevant towards a realization of 
such a proposal, without invoking the anthropic principle and without acausalities. 

To finish this discussion, we remark that in the presence of superluminal particles one 
can find examples where a stable system appears to have a tachyonic instability in some 
Lorentzian reference frame. It is straightforward to check that this happens, for instance, if 
one adds a mass term — l(m 2 n 2 ) to the action (29) (for v 2 > 1). 

To summarize, since the goal of the present paper is to find a connection between a 
violation of the NEC and the presence of local instabilities, we will adopt a very weak criterion 
of stability. We will assume that if the system (13) possesses superluminal excitations one 
should not conclude that it is unstable after observing (apparent) instabilities in some specific 
Lorentzian frame. Instead one should check whether it is possible to find some frame where 
the system (13) is manifestly stable. If this frame exists, it should be used to quantize the 
system thus defining the notion of causality 9 . We leave aside the issue whether an effective 
theory with superluminal excitations can arise from a consistent microscopic theory. 

To conclude let us stress that in the absence of superluminal excitations the above pecu¬ 
liarity does not arise and all Lorentzian frames are equally good for the analysis of stability. 

3.2 Arbitrary number of flavors 

Let us now discuss what the stability criterion is for a multiflavor system of general form (8) 
with I, J — 1,..,, N. We still focus, for the moment, on two-dimensional systems. We work 
in the limit discussed in section 2 in which the kinetic matrix Lu(p) is a quadratic form of 
momentum, 

Lu{p) — . (34) 

Dispersion relations for on-shell states are determined by the condition 

det Ljj(p ) = 0 . (35) 

This equation defines a set of rays in the (u, k) plane. As before, we can rewrite it in terms 
of the slope (velocity) v as defined by eq. (15), 


det L/j(z/, a) — 0 , (36) 

9 In the above discussion (and in the following sections) we have restricted to Lorentzian frames the 
search of coordinate where canonical quantization leads to stability. It is easy to see that this restriction 
is unimportant. As we explained, the choice of canonical quantization corresponds to the choice of time 
coordinate t. Assuming a good choice of equal time slicing exists, the invariant distance in that frame will 
in general be ds 2 = adt 2 + Cidtdx 1 + gijdx z dxi . It is clear that by a reparametrization of the form t = ct', 
x l = CjX' J + c l t’ the metric can be made Lorentzian without changing the time slicing. 
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where Lu(v) = Ljj(p)/k 2 and a = ±1 characterizes the direction of the spatial momentum 
k. Let us consider eq. (36) at k > 0 10 . This equation has 2N roots. In analogy with the one 
flavor case, upon quantization half of these roots {v+i} will correspond to creation operators 
for particles with momentum k and energy 


u = v +I \k\ . 


(37) 


The second half of the roots {v~i} corresponds to annihilation operators for particles with 
momentum —k and energy 

u = —v_i\k\ . (38) 

One necessary condition for stability is that all roots v±j be real. 

To understand the issue of ghosts let us understand how the two sets, {v+i} and 
are characterized. For this purpose, notice that (cfr. the analysis in Ref. [13]) at any value of 
z/, the real and symmetric matrix Ljj(v) has N real eigenvalues \i(v). When v takes one of 
the values in {v + i} or in {v-i} one of the eigenvalues A/ is equal to zero. We will now prove 
that the two sets of zeros are characterized by the sign of the derivative of the corresponding 
eigenvalue at its zero, 

d\j I 


and 


dv 

d\j 


> 0 if A iiy+j) = 0 , 


v +j 


dv 


< 0 if \t(v—t) — 0 . 


(39) 


(40) 


It is straightforward to check that these conditions agree with the results of subsection 3.1 in 
the case of one flavor. Conditions (39), (40) are easily derived by using the Kallen-Lehmann 
representation of the propagator. Consider indeed the Feynman correlator 

C"(() = (0|T( I r'(«) Ir J (0)) |0> 

and rewrite it by inserting a complete set of states | n) between n 1 (t) and vr J (0) 

G IJ < 




n 


—iE n t ji ,rIJ 
1 V 1 (n) 


E 

n 


9(—t)e 


,iE n t. i\,rJI 
1 V 1 {n) i 


where E n is the energy of the n-th state and 


= (0|^(0)|n)(n|7r J (0)|0) . 


(41) 


(42) 


Note that all the matrices are positive semi-definite, eq. (41) implies that the Fourier 
transform G IJ (E) of the propagator takes the form 

i 


G IJ (E ) = 




^ E- E n + ie M & 


^2 E + E n - ie MJ 2 


(43) 


10 Note that Lu( v,l) = 1), so that Lu(y, 1) alone carries all relevant information about the 

spectrum. 
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Consequently, the poles of —iG IJ (E) with positive residue correspond to the energies of the 
physical states. On the other hand, from the Lagrangian we have 

—iG IJ (E) = Ljj(E) . (44) 

Consequently, using matrix notations we may write 

-iG(E) = k 2 0 T (v)D(v)0(v) , (45) 

where D(y ) is a diagonal matrix with eigenvalues equal to Xj 1 ^) and 0(v) is an ortogonal 
matrix. Any pole of G(E) corresponds to a zero of some A 7 (z/), while the sign of the residue 
is the same as the sign of the first derivative d.X 1 (is)/dv. This implies eqs. (39, 40). 

Using these conditions one can determine the spectrum of the system. According to 
our previous discussion, the system is stable provided we can so define time that upon 
quantization there are no multiparticle states in the lower light cone (see eq. (26)). 

3.3 Arbitrary number of dimensions 

Let us finally discuss how the above results generalize to held theories in arbitrary (d + 1) 
dimensions. A straightforward generalization of the dispersion relation (36) is 

det Ljj(u, n) = 0 , (46) 

where, as before, Ljj(u, n) = Ljj(p)/ k 2 and n is now a unit d-dimensional vector in the 
direction of the spatial momentum k. 

Again the system (8) is stable if there is a Lorentzian reference frame, where eq. (46) 
has 2 N real roots for any n and there are no states in the lower light cone. Note that our 
derivation of the conditions (39), (40) was not sensitive to the number of the dimensions. 
Consequently, these conditions are applicable for arbitrary d as well. The only difference is 
that now roots are functions of the unit vector n, and iy + (n)-roots correspond to the 
particles with momentum along n while ///_ (n)-roots correspond to particles with momentum 
in the opposite direction. 

4 NEC violation and instabilities 

To understand the implications of the NEC violation for stability, let us first study how 
stability itself is encoded in the Lagrangian. We start with the analysis of two-dimensional 
systems. Consider a stable system of the type (1). According to the discussion of section 3 
this means that a Lorentzian reference frame exists where the spectrum of states does not 
contain either tachyons or states with energy-momentum in the lower light cone. Let us 
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see in more detail how this restricts the properties of the matrix Ljj(p) in eq. (8). As we 
discussed in subsection 3.2 the properties of the spectrum can conveniently be read off from 
the properties of the eigenvalues A i(u) of the matrix 

Lu{u) = L u (p)/k 2 , (47) 

where without loss of generality we assumed that k > 0. Absence of tachyons implies that 
these eigenvalues have 2N zeros in total. Let us assume that we are in a generic situation, 
where all these roots are simple. Multiple roots may emerge due to some coincidental 
degeneracy or as the result of a symmetry. We postpone the discussion of degeneracies 
to section 6. 

Let us now understand what root configurations on the u -axis are implied by stability. 
As before, we indicate by {z/ + } and {z/_} the sets of roots with respectively positive and 
negative slope. In what follows we consider all the particles together, so that we include in 
the sets {v±\ the roots v = ±1 associated to the “graviton”, “photon”, etc. By eqs. (37 - 
40), if there exist zq + G {zq + } and zq_ G {zq_} such that Ui + < then we can construct 
a two-particle state with zero total momentum and negative energy 

E = (u 1+ - v x _)k . (48) 

By the discussion in the previous section we thus have a ghost instability. On the other hand 
it is straightforward to check that if 

vt+ > v.j- (49) 

for any two members of the two sets (see Fig. 5), then the Hilbert space does not contain 
any state with zero energy-momentum other than the vacuum, so that there is no ghost 
instability, eq. (49) is just a generalization of the conditions (37), (38) obtained in the one 
flavor case. 

Notice that eq. (49) implies that no eigenvalue A i(v) has more than two zeroes on the v 
line. Since we assume that there are 2 A" zeroes in total (we barr degeneracies), then each 
eigenvalue must have precisely two zeroes. The behaviour of the Athat follows is shown 
in Fig. 5. Since the signs of the eigenvalues define the sign properties of the matrix Ljj(u), 
we find that stability holds if and only if the following two conditions hold: 

(i) There exists z/ max > 1, such that Ljjiy) is positive definite for \u\ > // max • 

(ii) There exists —1 < v m \ n < 1, such that Ljj(u min ) is negative definite. 

Condition (i) is just the usual requirement that the kinetic part of the Lagrangian (8) be 
positive definite. Moreover if condition (ii) is satisfied for the special value z/ m j n = 0, then 
the purely gradient part of the Lagrangian (8) is negative definite. Conditions (i), (ii) with 
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zv ni in = 0 are equivalent to positivity of the Hamiltonian. In general, however, validity of (i) 
and (ii) is a weaker request than positivity of the Hamiltonian. 

At this point it is convenient to go back to the matrix Ljj(p ) and reformulate conditions 

(i) and (ii) in a manifestly Lorentz covariant way. Let us define the following scalar function 
of 2-vectors p and q 

L(p , q) = Lu{p) A^f A J u q u . (50) 

In practice we are evaluating the quadratic Lagrangian of a perturbation 

V = . (51) 

Perturbations of this special form can be obtained by acting on the vacuum (7) with a 
space-time translation + q /J . Perturbations (51) are therefore just the Goldstones 

(sound waves) of broken space-time translations. That is why we do not lose any relevant 
information by focusing on these perturbations. By using the defining eqs. (9, 11), we obtain 
a simple relation 

L(p, q) - L(q, p) = \ (p 2 T^q v - q 2 T^p u ) = f(q, p, T, w ) , (52) 

which relates the kinetic matrix to the energy momentum tensor, and which will prove crucial 
in relating NEC to stability. As we shall see shortly, stability in the form of conditions (i), 

(ii) forces the function / to be positive for properly chosen p M and q A On the other hand, 
if NEC is violated then T^ v is negative definite along some direction. Our strategy will be 
to prove that under generic conditions there exists a tension between these two facts. 

Now, notice that \v\ > 1 and \v\ < 1 correspond respectively to time-like and space-like 
2-momentum. Conditions (i), (ii), then imply the following two properties of L(p,q). First, 



Figure 5: The eigenvalues A j of Ljj as functions of u. The configuration shown corresponds 
to a stable system: all the positive derivative roots ui + are on the right of the negative 
derivative roots (taking into account also the roots ±1 corresponding to the dispersion 
relation of the graviton). As clear from the picture, this property is equivalent to conditions 
(i), (ii). 
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there exists a cone T + in the region of time-like momenta which includes the (u 7 ^ 0 , k = 0 )- 
axis such that for any f eT + and for arbitrary q 

L(t,q)> 0. (53) 

Second, there exists a space-like line T_ such that for G T_ 

L(s,q)< 0 (54) 

for any q^. By eqs. (53, 54) it then follows that 

f{s,t,T liV )> 0, (55) 


for t G T + and s G T_. 

It is now straightforward to extend the discussion to systems in d > 2. Consider indeed 
the subspace of the states where all particles move collinearly (he the spatial momentum is 
along some given direction n). Then one can literally repeat the arguments above and obtain 
that the conditions (i), (ii) should hold for any direction n. (In analogy with subsection 
3.3 the matrix Lu{y) now depends on n as well). We stress however that these conditions, 
while necessary, are no longer sufficient to ensure the absence of ghost instabilities. Indeed, 
for anisotropic systems one may imagine a situation where states with total momentum in 
the lower light cone are absent if one considers only states with collinear particles, but are 
present in general. Nonetheless, these necessary conditions allow us to establish a relation 
between NEC and stability in a rather broad class of theories. On the other hand we will 
be able to construct examples of stable systems that violate NEC. Therefore we do not see 
much motivation in trying to formulate stronger conditions. 

Focusing just on states of collinearly moving particles the conditions (i), (ii) are easily 
generalized. As a consequence condition (53) is unchanged: there must exist a cone T + in 
momentum space, which contains the cn-axis and such that the matrix Ljj(p ) is positive 
def ini te for any p G T + . Condition (54) is similarly generalized: as before, for any given 
space direction n there exists a line of space-like momenta along which the matrix Ljj(p) 
is negative definite. By varying n we find a skirt-shaped cone T_ (see fig. 6). As before 
stability implies that eq. (55) holds for t G T + and s G T_. 

Now that we have characterized the implications of stability, let us assume that the NEC 
condition is violated along some null vector rffi: T^ v <0. In d > 2 by continuity 
there must exists a whole neighborhood of null directions along which NEC is violated. The 
two-dimensional case is singular in this respect, since there the light-cone is made up of two 
disconnected directions; we will analyze the implications of this fact later. Let us therefore 
concentrate on d > 2 for the moment. Any pair of such linearly independent NEC-violating 
null vectors n 1 and 77,2 defines a plane A f. Now, both the NEC and the scalar quantity / are 
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Figure 6: The various regions in the (cu, k ) space discussed in the text. Ljj is positive dehnite 
inside the yellow narrow cone T + , while it is negative dehnite on the skirt-shaped blue cone 
r_. The green plane J\f is one of the planes on which the stress-energy tensor can be made 
negative dehnite. The grid-like cones are the light-cones. 

invariant under the shift symmetry 

T/ja, -> = T ta , + a rjftv , (56) 

for arbitrary a. By choosing a = —T ]W n^n 2 /(n\ ■ 77 - 2 ) we can make strictly negative 
dehnite on the whole plane Af, i.e., on all linear combinations of n 1 and n^- Notice that 
M is a time-like plane: it intersects the light-cone along two lines (defined by ri\ and n 2 
themselves), while a space-like plane would only intersect the light-cone at the origin. By 
varying ri\ and n 2 inside the ‘NEC-violating neighborhood’ we can construct a whole family 
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of such time-like planes A f on which T^ u can be made strictly negative by means of the 
suitable plane-dependent shift. Then if we can prove that at least one of these A/"’s must 
intersect both T + and T_ we have a contradiction. Indeed, on the one hand stability forces 
/(s, £, T M „) to be positive for all s e T_ and t e T + . On the other hand, if NEC is violated 
and one of the A/"’ s intersects both T + and T_ we can pick t M and s M along these intersections. 
In this case T /IU is negative definite both on F and on and /(s,t,T Ml/ ) is thus negative, 
by its very definition eq. (55). 

It is evident that all A/"’s must intersect T_, since the former are time-like planes while 
the latter is a space-like hypersurface. On the other hand they may not intersect T + . This 
clearly happens when, in the frame where stability is manifest, all A/"’s are almost tangent 
to the light cone (which is the case if the neighborhood of NEC-violating null vectors is 
very small) and T + is a tiny cone around the energy axis. This suggests that in order to 
establish a connection between NEC and stability we need extra assumptions. These are 
indeed offered in at least two physically relevant situations: 

• The first is the case where all excitations are subluminal, for which the region T + 
necessarily contains the light cone. Then all planes J\f intersects T + and we obtain 
that the assumption of stability and NEC violation contradict each other. 

• The other example consists of systems in which superluminal modes may be present 
but there are two NEC violating null vectors with opposite spatial components rif 2 = 
(l,±n). In this case the plane A f generated by ri\ and n 2 contains the cu-axis and 
necessarily intersect the T + -region. This property trivially holds for systems, where 
the NEC is violated for all null vectors, in particular for isotropic systems. This is the 
case which is relevant for cosmological applications. 

The fact that NEC follows from stability and subluminality, is perhaps not completely 
surprising. As discussed for instance if ref. [20] the dominant energy condition, which implies 
NEC, corresponds to a background with positive energy density and subluminal energy flow. 
However the energy density and flow of the background do not precisely correspond to the 
Hamiltonian and velocity of the fluctuations around the background itself. Mathematically 
this is seen by comparing eq. (9) and eq. (11). The matrix describing the propagation of 
perturbations involves also the second derivative tensor Fjjkl while the energy momentum 
does not. Conversely T^ v involves the cosmological constant term oc g tJiV F which does not 
affect the propagation of the modes. These differences suggest that a quick conclusion that 
stability and subluminality of the modes implies NEC would be too naive. However, the 
quantities we have used in eq. (52) somewhat cleverly eliminate the contributions of F IJKL 
and F in such a way that the naive conclusion is correct. 

Let us now come back to the two-dimensional case. In d = 2 the light-cone is composed 
of the two null directions n± = (1, ±1). If NEC is violated along both null directions, then, 
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as before, the stress-energy tensor can be made negative definite on the plane J\f generated 
by n .|_ and n_, which in this case coincides with the whole (u, k ) plane! This plane obviously 
intersects (in fact contains) both T + and T_. We are therefore led to the above contradiction. 
If instead NEC is only violated along, say, n+, then the region A/" on which the stress-energy 
tensor can be made negative definite is just the line generated by n+. In this case J\f cannot 
intersect both T + and T_, and there is no evident contradiction. A peculiar feature of two- 
dimensional systems is that the condition that NEC be violated along one and only one 
null direction is equivalent to the requirement that the energy-momentum tensor cannot be 
diagonalized by a Lorentz-boost. Indeed they are both equivalent to the inequality 

iToo + Tnl < 2|T 0 i| . (57) 

We can therefore summarize the d = 2 case as follows. If the stress-energy tensor is cliago- 
nalizable by a Lorentz boost and the NEC condition is violated, then the system is unstable. 
On the other hand we will now show that by giving up the assumption that 7j tJ , be diago- 
nalizable we can construct a system that violates NEC while being perfectly stable, indeed 
with a positive Hamiltonian. We stress again that such an assumption is only crucial in 
d = 2, given its accidental relation to the ‘extent’ at which the NEC is violated. Indeed in 
the next subsection we will also show an example in d > 2 of a stable NEC violating system 
with diagonalizable energy-momentum tensor. 


4.1 Examples of stable, NEC-violating systems 

In c/ = 2, consider two scalar fields (j ) 1 , / = 1,2 with a purely derivative action and choose 
the following classical background 


4>° = t, (f) 1 = x. 


(58) 


By defining 

B IJ = B IJ - 5 IJ (59) 

(with the “meson” held B IJ defined in eq. (2)) we now present an action written in term 
of B 1J such that T jiv satisfies eq. (57) and such that the quadratic Hamiltonian of small 
perturbations is positive definite. The advantage of the definition (59) is that B IJ = 0 on 
the background (58). This simplifies some of the algebra. The Lagrangian is of the form 


£ = (H 00 + B 11 + 2(1 + a 2 ) B 01 - (a 2 + b 2 ) B m B m + (a 2 + c 2 ) B u B 01 ^j . (60) 


with 

b 2 < 2, c 2 < 2 . 
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The energy-momentum tensor for this system can be calculated using equation (11), 


T 

± /J,1S 


( 1 1 + a 2 \ 

\1 + a 2 1 ) 


( 61 ) 


which satisfies (57) and violates the NEC. To study the stability of the system we need to 
write the Lagrangian for the fluctuations and this is done in the appendix. 

To conclude this example, let us note that this construction can be readily extended to 
provide an example of a stable system with non-diagonalizable (and, consequently, violating 
NEC) energy-momentum tensor for an arbitrary number of dimensions d. In addition to 0°, 
0 1 one introduces more scalar fields 0*, i — 2 r ..., d — 1 with background values 


til = x* • 


As before we will denote small perturbations of the 0*-fields about this background by 7r*. 
We also add the following term to the Lagrangian (60) 

AC = (n 2 B h - f 2 [( B 0i ) 2 + ( B li ) 2 ]) (62) 

where we generalized the definition of the B{j in the obvious way. For large enough f 2 the 
last two terms in eq. (62) provide a dominant contribution of the correct sign to the gradient 
terms of the 7r 0,1 -fields in the ad-direction. The first term in eq. (62) provides healthy kinetic 
and gradient terms for the 7T*-per turbations. In the limit N 2 3> a 2 , b 2 , c 2 , 1 one may neglect 
extra mixings involving 7r'-fields and coming from the last two terms. 

As anticipated, we now want to show that the requirement that T /iiy cannot be diagonal¬ 
ized by a Lorentz boost is not crucial at d > 2. Indeed we construct an example of a stable 
NEC violating system with diagonalizablc energy-momentum tensor. We work in d — 3, the 
generalization to d > 3 being straightforward. As a first step, we take the stress tensor, in 
the frame where it is diagonal, to be of the form 

/a 2 0 0 \ 

- 0 -b 2 0 (63) 

V o o l / 


with (a 2 — b 2 ) < 0. By the discussion in the previous section the reference frame where a 
system with such X0, has a chance of being manifestly stable, must be boosted in such a way 
that all the null NEC-violating vectors are close to each other and the jV-planes they define 
are thus almost tangent to the light cone. The above stress-energy tensor violates NEC for 
null vectors whose spatial component mainly extends in the x 1 -direction. We thus consider a 
large boost in the x 2 -direction. In the boosted frame the energy-momentum tensor becomes 

/ ch 2 a 2 + sh 2 0 chsh(l + a 2 ) \ 

T^=\ 0 —b 2 0 (64) 

\ ch sh (1 + a 2 ) 0 ch 2 + sh 2 a 2 ) 
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where sh = sinhA and ch = cosh A 3> a 2 , b 2 , 1, are the boost parameters. To construct a 
stable system with stress tensor given by eq. (64), we introduce three scalar held (ft 1 with 
vev’s 

€ = x 1 . (65) 

and write the Lagrangian in term of the mesons B 1J defined previously 

C = T U B IJ + (a 2 [(B 01 ) 2 - ( B n ) 2 - ( B 12 ) 2 ] + 

chsh (1 + a 2 ) [B 02 (-B 11 + B 22 ) + 2R 01 R 12 ]) (66) 

with a 2 ~ b 2 ~ a 2 ~ 1. The stability of this system is discussed in the appendix. In 
agreement with our general expectations it contains a highly superluminal mode. 

5 Effective field theory for solids and fluids 

In the previous sections we focused our attention on systems involving derivatively coupled 
scalar fields. It is not obvious a priori whether our results can be exported to different physical 
systems. In particular, can we apply them to fluids? After all, most of the cosmic history is 
dominated by a mixture of fluids, their physical relevance in cosmology certainly overwhelms 
that of scalar fields. In this section we show that our results are straightforwardly applied 
to universe filling continua, solids, and in particular to perfect fluids. Indeed the dynamics 
of these systems is described by an effective theory of three derivatively coupled scalars, as 
we will now discuss. Although most of the results we will discuss are well known, we believe 
our presentation to be still useful as it provides a direct connection to the subject of the 
previous sections and to the investigations on modifications of gravity. For a recent thorough 
review on the Lagrangian formulation of fluid dynamics we refer to Ref. [21], Our treatment 
is closer to the effective (non-relativistic) theory for supersolids developed in Ref. [22], but 
we work from the beginning with the fully relativistic formulation. 

5.1 Generalities 

Let us consider a mechanical continuum filling all of (3-dimensional) flat space and such 
that at each space point there resides just one (and only one) dynamical entity. There 
thus exists an invertible mapping <ft between coordinate space and the internal space A4 
labeling the dynamical variables (for instance the positions of the atoms in a crystal). This 
is made explicit by working with a held theory of three scalar fields <ft J (x, t), I = 1,2,3, 
which describe both the mapping and the dynamics. Since A4 is diffeomorphic to M 3 we 
can always choose held coordinates (ft 1 such that in some given stationary state, for instance 
the ground state itself, the mapping is just 0 7 (aA,f) = x 1 . In what follows we shall indeed 
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assume 0 7 (a7,f) = x 1 on the ground state. Notice also that the 0 7, s correspond to a choice 
of space coordinates comoving with the solid. This is to say that, given an arbitrary time 
dependent configuration 0 7 (V,t), we can choose space-time coordinates by 

x n = 4> l (x, t) t' — t (67) 

for which the scalar fields take constant value 

(j)' 1 {x', t') = 4> I (x(x'), t') = x' 1 . (68) 

Of course in the primed coordinates all the information about the state of the solid is 
contained in the metric, which is no longer Minkowskian. 

Let us work in Minkowskian coordinates around the ground state 0 7 (ad,f) = x 1 , and 
let us consider the symmetries of the system. Ordinary space translations T x : t) —> 

(j) 1 (x + a, t) are obviously a (spontaneously) broken symmetry. On the other hand we would 
like to focus on homogeneous solids, for which there still exists an unbroken translational 
invariance. It is evident that the only way to achieve that, is to further assume the additional 
internal shift symmetry 

T^\ ¥{ x i t) —> ^{x, t) + c 1 . (69) 

This way the ground state is invariant under a combination of ordinary translations and 
with a 1 = — c 1 . Therefore T x — represents the conserved three momentum for processes 
corresponding to small fluctuations around the vacuum (for instance in the scattering of 
phonons). The broken generators T x + T ^ are instead in one to one correspondence, by 
Goldstone theorem, with 3 gapless excitations, the phonons. At the lowest derivative level, 
the relativistic invariant Lagrangian for such a system will clearly be written as a function 
of the ‘meson’ fields B 1J = d^(j) 1 (j) J . It thus falls precisely into the class of held theories 
we have been discussing throughout the paper! If the Lagrangian is a generic function of 
the B 1J , the system will describe an homogeneous but otherwise anisotropic solid. More 
symmetric systems can be constructed by assuming additional internal symmetries acting 
on (p 1 . For instance we can consider an invariance 

¥ - D 7 j0 J (70) 

where D 1 j belong to some discrete subgroup of 0(3). In this case our system will have the 
same macroscopic properties as a crystal with the corresponding symmetry. We can also 
consider the limiting case in which D 1 j is extended to the full rotation group SO( 3). In 
which case our system will essentially describe a jelly. 

The jelly system does not correspond yet to the most symmetric possibility. Indeed one 
can consider a system for which this symmetry is extended to the group of volume preserving 
diffeomorphisms 

¥ -»■ 0 / (</>), det = 1 . (71) 
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While this is a huge internal symmetry it still allows a non trivial invariant quantity: the 
determinant of B IJ . The most general Lagrangian for such a system is then 

C = F(B) , B = det B IJ , (72) 

where F is a generic function. It is pretty obvious that, with three fields, we cannot further 
extend the volume preserving diffeomorphism to a larger symmetry that still allows for a 
non-trivial action. What is the physical interpretation of a system of this type? The answer is 
intuitively obvious: our Lagrangian describes the dynamics of a (relativistic) fluid. Consider 
indeed the ground state of the system 0 J (x) = x 1 . We can imagine, starting from this 
configuration, to build another one by moving around the individual points: x —> f(x). In 
the new configuration the field profile will be given by 

= (r 1 ) 7 (x) • (73) 

Now, the intuitive notion of fluid suggests that its energy does not change if we do not 
compress it, no matter what different form its individual parts assume. Therefore we expect 
that for a coordinate transformation / that preserves the volume in physical space, the 
energy of the system should remain the same. This is obviously true in our case because 
eq. (73) can be compensated by an internal symmetry transformation. 

5.2 Fluids 

In this subsection we want to work out in some detail the properties of the relativistic fluid 
that follow from the Lagrangian in eq. (72). We want first to make contact with the standard 
description that can be found, for instance, in Ref. [23]. In particular we want to find the 
field theoretic description of quantities like the fluid 4-velocity, the current, the pressure and 
energy density. 

The definition of the 0’s as comoving coordinates makes the characterization of the fluid 
velocity field u^(x) immediate: is a vector field along which all the 0 7, s stay constant, 

(x(t)) = uFd^ 1 = 0 1 = 1, 2,3. (74) 

dr 

These orthogonality constraints together with the normalization constraint uJ'u^ = —1 com¬ 
pletely determine u M (except for the discrete choice of the direction, which is set by the 
request that be future-directed). Indeed the orthogonality constraints forces to be 
proportional to the ‘vector product’ e ^ 7 d^dg^d ^ 3 oc e MQ/37 e IJK d^<\) K . The 

proper normalization is 

^ = -^7= ^ e TJK . (75) 
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Note that this expression is invariant under generic (not necessarily volume-preserving) dif- 
feomorphisms acting on the 0’s. This is correct, since what matters for characterizing is 
only that the 0’s are comoving coordinates. By the above equation, is associated to the 
volume form 12 over the fluid space 


u l jdx 1 ' 1 


1 

71 


★ (\euKdct) 1 A d(j) J A d(j ) K ) 


4=*n 

Vb 


(76) 


and the one form *12 = J^dx 11 is naturally interpreted as the current of fluid points. Consis¬ 
tent with this interpretation, the current is conserved 


d -k (J^dx^) = d^ejjKdcp 1 A d(j) J A d(f> K ) = 0 . 


(77) 


Notice that the conservation of J /t , rather than arising on-shell from a symmetry of the 
Lagrangian via Noether’s theorem, is already valid off-shell since the number of degrees of 
freedom, that is the volume of A4, is fixed by construction. From J^ — n- it also follows 
that the density of particles n is 

n = Vb (78) 

What about the equation of state? The stress-energy tensor is given by eq. (11), with 

A/J - ~qbTT ~ A U t> TJ , 

T IW = 2 F' B By} - g^uF . (79) 


On the other hand for a perfect fluid the general expression is = (p + p) + pg^ u . 
The energy density and pressure therefore can be extracted by contracting with g^ v and 
with u^u v . By making use of the constraints eq. (74) we get p and p as functions of B, 


p — —F{B) , p = ~2F'(B) B + F(B) . 


(80) 


Note that although the stress-energy tensor eq. (79) for a generic configuration depends 
on the full structure of d ^ 1 , the energy density and pressure depend on B only. This 
is physically correct: like the Lagrangian they must be scalar quantities only sensitive to 
compressional modes. 

By eqs. (78, 80), the three quantities n, p, p are determined by just one parameter. On the 
other hand, the equation of state of a fluid corresponds to just one constraint f(n,p, p) = 0. 
The additional local constraint arises, in our description, from the homogeneity of the initial 
entropy density together with the isoentropicity of the fluid transformations, as we now 
explain. In the ordinary treatment, see [23], the current and energy momentum conservation 
equations (<9 /i J /i = d lL T tLV = 0) can be combined to give an equation that establishes the 
conservation of the entropy per comoving volume a along the flux 

0 = vfdpo = -^{pd^l/n) + d^p/n)}^ (81) 
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In other words the motion of each fluid portion is isoentropic, although the entropy density 
a varies in general through space. Therefore the above equation introduces an additional, 
non-algebraic—he., differential—constraint between n,p,p. Now, the fluid we are describing 
by our held theory is homogeneous from all stand points, in the sense that any given con¬ 
figuration can always be brought to the homogeneous ground state by means of mechanical 
transformations. In other words, in our formalism there is no room for genuine thermo¬ 
dynamic degrees of freedom: to increase the temperature we have to compress the fluid, 
we cannot ‘heat it up’. In particular our fluid is homoentropic (the entropy per comoving 
volume is the same for all fluid elements), in which case the last equation is replaced by the 
stronger one <9 /( cr = 0 which corresponds precisely to the additional algebraic constraint we 
have between n,p, p. 

We can specify our analysis of quadratic fluctuation to the fluid. Consider a background 
(j) 1 = ax 1 1 where a parameterizes the level of compression B = a 6 , and let us parameterize 
the perturbations by (j) 1 = a{x : + it 1 ). By expanding eq. (72), integrating by parts and 
neglecting total derivatives we get the Lagrangian for the Goldstones, 

C = —F'(B)B ■ (tt 7 ) 2 + (F'(B)B + 2 F"(B)B 2 ) • (d/Tf 7 ) 2 . (82) 

Only the longitudinal Goldstone has a gradient energy term; as expected it is the only propa¬ 
gating degree of freedom: sound waves are longitudinal. The dispersion relation corresponds 
to a speed of sound c| 

F’ + 2F"BB p’(B) dp 
° s ~ ¥' = p’(B) ~ Tp 

where in the last equality we have made use of the equations of state eq. (80) and have thus 
made contact with the usual expression for the speed of adiabatic sound waves in a fluid. 
On the other hand, the transverse phonons, 7 f oc e l ^ k djfk do not propagate. Their general 
solution at the infinitesimal level is linear in time 

Ti l r = e^dj (, f k (x ) + g k (x)t) (84) 

and corresponds to vortices in constant rotation. In the linearized approximation the fluid 4- 
velocity is indeed id 1 = (1, — 7r*), so that we have 11 v l = e^ k djg k {x). The trivial time evolution 
of infinitesimal vortex perturbations follows directly from the internal volume preserving 
diffeomorphism invariance of the fluid, eq. (71). As a matter of fact, (also at the full non¬ 
linear level) this symmetry implies the existence of an infinite set of conserved currents and 

n For large t the linearized approximation for this solution breaks down. Notice however that the bulk of 
the value of tt’ t at large t can be eliminated by a volume preserving diffeomorphism around the stationary 
solution. This means that what matters physically is just the time derivative t: 1 t , which is just the fluid local 
velocity. 


(83) 

<r=c,onst 
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of the associated charges 


= (85) 

where f 1 is any vector function satisfying df 1 /dtp 1 = 0. The conservation of the infinite set 
of associated charges basically corresponds to the trivial dynamical evolution of the vortices. 

This can be more explicitly seen by working in the comoving coordinate frame, i.e. by 
choosing the p 1 as space coordinates and by taking x l ((ft,t) as the dynamical variable. From 
the basic relation (ft 1 (x((ft,t),t) = const we obtain by differentiation 

ftft 1 = -dktfx* = -W\ x k . (86) 

By noticing that det^*- 7 — id id) = (1 — x 2 ) we have 

B = (detkF) 2 (l - x 2 ) , (87) 


so that the action now reads 

S 


dtd 3 (ft 


F ((detkF) 2 (l 
|detW| 



( 88 ) 


This action is manifestly invariant under volume preserving diffeomorphisms acting on the 
comoving coordinates (ft 1 . On the under hand Lorentz symmetry acts now as a combined 
transformation of the time variable t and the dynamical fields x % ((ft,t), and invariance of the 
action is not explicitly manifest. The volume preserving cliffs are realized on the dynamical 
variables as 8x l = f I dix l , with djf 1 = 0, so that the set of conserved charges is given by 

C f = I d s pf- f T d T F (89) 


For f 1 = a 1 = const the 3 conserved charges just correspond to spatial momentum. The 
genuine fluid invariants are obtained by choosing f 1 = e IJK dj(a,K8 3 ((ft — 0o))? for arbitrary 
constant ax and 0 O . By eqs. (88, 89) we obtain, after integration by parts, that the local 
2-form 


d. 


F^detlF) Xidxx 1 d<ft J A d(ft 


I< 


d{F'VB Ui dx*) 


(90) 


is time independent, where in the last equality we have used the relation x t = \/l — x 2 u t for 
the space components of the 4-velocity. Considering any 2-surface E bounded by a closed 
circuit C = <9E we have then 


-f- (£ F'VB Uidx 1 = — [ d 
dt Jq dt J s 


F'VB Uidx' 1 


= 0 


(91) 


From the point of view of coordinate space the circuit C depends on time as it moves with 
the fluid flux, so the statement is that the line integral of F'y/Bui is conserved along the 



flux. This quantity computed for all closed paths C measures indeed the vorticity of the 
fluid’s configuration. Onr result states that vorticity, thus defined, does not evolve with 
time. This result corresponds to the relativistic generalization of Kelvin’s theorem. In the 
non-relativistic limit, we have p <C p, which from eq. (80) corresponds to F(B) ~ \/~B. In 
this case the conserved quantity becomes simply the circuit integral of the velocity as 
stated in Kelvin’s theorem. 

In the last discussion Lorentz invariance was not explicitly manifest. We can however 
easily express the result in eq. (91) in a manifestly Lorentz covariant manner. From the 
definition of the 4-velocity and from eq. (91) the relevant one-form we should focus on is 

F'VB u^dx' 1 = V tl dx 11 . (92) 

Using V the equations of motion 

d fi [F'D Bjj d' l (P J ] = 0 (93) 

are written conveniently as 

(dV) A dc[) J A d(j) K e IJK = 0 (94) 

The geometric implication of this formula is that (dV) = 0 on any two surface E tangent at 
all points to the four velocity This is easily seen by working in comoving coordinates, 
where u M = (1, 0, 0, 0) and eq. (94) reads (dV) 0I = 0 for / = 1, 2, 3. We can create such 
a surface E by choosing a spacclikc closed cycle C 0 and by considering the world-surface it 
describes along the fluid motion up to some final cycle C\. We have <9E = Co — C\ so that 
by Stokes theorem we find 

0 = (f dV = (f V - l V (95) 

is Jco JC\ 

As expected, the quantity §F'\Z~B u^dx^ 1 is conserved along the flow. Notice however that 
neither Co nor C\ has to be necessarily chosen at fixed time. To summarize, our result states 
that given a two-surface E with the topology of a cylinder and tangent to the local 4-velocity 
at any point, the integral of V = F'\/~B u^dx^ 1 around E is independent of the cycle chosen 
to compute the integral. 

5.3 From Fluid to Superfluid and to Ghost Condensate 

We have seen that vorticity is conserved through the time evolution of the system, corre¬ 
sponding to the ‘trivial’ dynamics of the vortices. In particular if vorticity vanishes on any 
initial time surface it will vanish throughout the time evolution. For this subclass of the 
fluids configuration it is natural to expect our system to be described by just one scalar held 
instead of three. Let us see how this works. 
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The most direct, although at first sight unintuitive, way to proceed is to introduce an 
auxiliary 1-form V0 and start with the following action 


S vfl = / G{X) — V A 


(96) 


where, X = —V^V^ and hi is the fluid volume form defined before. The V equations of 
motion 

-2 G'{X)V» = (★ft)'* = J' 1 (97) 

allows us to algebraically solve for V) t as a function of and to write eq. (96) as a functional 
of the fluid variables 


G(X) -V AQ 


= / F(B). 


(98) 


I SS/SV =0 

The functions G(X = — V 2 ) and F(B = — J 2 ) are the Legendre transforms of one another, 
and Vj, and are the conjugated variables. By the usual properties of Legendre transforms, 
eqs. (96, 98) also imply the conjugate to eq. (97) 


2 F , (B)J li = V li . 


(99) 


Up to a factor of 2, the form V coincides on-shell with the form V we defined in the previous 
section, hence our notation. (This is also evident by varying eq. (96) with respect to 0/ and 
comparing to eq. (94)). Notice that eqs. (97, 99) imply the scalar relation 


2 F\B) = 


as well as 


2 G'(X) 

X = B(F\B)f 


B = X(G\X)) 2 
Moreover eq. (98) can be also written as 

G-F = -2 BF' = 2 XG'. 


( 100 ) 

( 101 ) 

( 102 ) 


We have now all the ingredients to discuss fluid motion at zero vorticity. From the 
discussion at the end of last section, this situation is realized when 


V = dy (103) 

where 0 is a scalar field. In this case oc C00, corresponding to the Lorentz covariant 
generalization of potential flow in a fluid [24], When eq. (103) holds, the equation of motion 
for 0 are trivially satisfied. However the dynamics of 0 is non-trivial and is determined by 
taking the divergence of eq. (97) 

d^G'(X)d^) = 0. (104) 
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We have therefore established that the solutions in the absence of vorticity are into one- 
to-one correspondence with the solutions for the dynamics of a derivatively coupled scalar 
with Lagrangian L = Moreover since d^ oc is a time-like directed vector, 

the scalar ip depends linearly on time in the ground state ip = ct. This is precisely the 
dependence of the charged held phase in a superfluid. We have thus established that, in the 
absence of vortices, the classical dynamics of a fluid is equivalent to that of a superfluid with 
suitably chosen Lagrangian (cf. [25]). As a further check of this equivalence, notice that the 
energy momentum tensors in the two descriptions coincide. This is manifest by eq. (98), 
since the fluid and superfluid action differ by the integral of the four form V Ail which does 
not depend on the metric. It can also be checked directly by using eq. (102) 12 . 

The relation between fluid and superfluid is a close relative of the duality between a 
2- and a 0-form. Consider indeed eq. (96), with hi = dA 2 . The 2-form A 2 acts like a 
Lagrange multiplier setting dV = 0 thus reducing to a 0-form V = dip action. On the 
other hand, by solving the V equations first we would get a fully equivalent local action for 
Ai- The fluid is distinguished from the two form because hi = dtp 1 A d<p 2 A d(p 3 is not the 
most general held strength a two form can have (it is a function of 3 held variables, instead 
of 6). However although the off-shell held space of the fluid is of smaller dimensionality 
that for the 2-form theory, the space of solutions is nonetheless bigger. This is perhaps 
counterintuitive, but it becames less so if one recalls that in a gauge theory some of the 
held variables act just as Lagrange multipliers. We do not want to explore these issues 
further here, but just make a few extra remarks elucidating the duality character of the 
huid/superhuid relation. Notice that dynamical equations and structural constraints are 
interchanged. So the conservation of J,j is structural for the fluid while it follows from the 
equations of motion in the supehuid. Conversely vorticity is dynamically conserved in the 
fluid and structurally zero in the superfluid. More importantly the fluid/superfluid relation 
is only valid for the classical action on-shell, and breaks down when quantum corrections are 
taken into account. In fact in a duality relationship the quantum perturbation series are not 
mapped into one another, and when one description is weakly coupled that other is not. 

This last property and others are made clear by focussing on a particular superfluid, 
the so called ghost-condensate model, which also provides an interesting and consistent 
modification of gravity [11]. The model is defined by assuming that the function G is 
positive and has a minimum (where C = 0) at some finite point X = c 2 in held space. This 
is the ghost-condensate point. At this point, although the held is time dependent ip = c*t, 
the energy momentum tensor is still proportional to the metric. Thus the universe is in a 
de Sitter phase with an effective cosmological constant = —G(c 2 ), which differs from the 

12 Indeed we originally derived our relation between fluid and supefluid by equating their energy-momentum 
and number current. This leads precisely to the relation described at the beginning of this section. We believe 
however the Lagrangian approach we followed is analitycally more transparent. 
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genuine cosmological constant defined at zero field —G(0). What does this point correspond 
to in the fluid description? By the first of eqs. (101) the answer is B = 0, i.e the point of 
infinite dilution, while, not surprisingly, eq. (100) tells us that this is a singular point. More 
precisely, by solving the second of eqs. (101) we find two solutions, corresponding to a branch 
singularity 

F ± (B) = -^/Bf ± (B) + F 0 (105) 

where f± are smooth functions satisfying 

lim f±(B) = ±c*. (106) 

B^O 

F + and F_ respectively corresponds to the superfluid branches X > c* ( G' > 0) and A" < cl 
(i G' < 0). As B —> 0 the F + branch becomes a diluted dust with positive energy density 
p ~ c*\fB and positive pressure 13 superposed to a cosmological constant F + (0) = F 0 . Notice 
that by eq. (102) F 0 = G(c*), which propels de Sitter expansion, is now truly the cosmological 
constant. On the ther hand, the branch F_ corresponds to a diluted fluid with negative 
energy density p ~ —c*\/£? but positive pressure p = B 3 / 2 f'_(B). The two branches X > cl 
and A" < cl thus correspond to two distinct fluids. These two fluids cannot mix as they 
are always separated by a region of infinite dilution, where B — 0. Moreover the regions 
corresponding to F + will tend to expand as they correspond to regions of higher density and 
pressure in an ordinary fluid. The regions corresponding to iA will instead contract [26, 27] 
as they correspond to a fluid with negative energy density but whose pressure grows with 
compression 14 . Notice also that in the region with G' < 0 the superfluid theory features 
a gradient instability and this is consistent with the imaginary speed of sound implied by 
eq. (82). However in the fluid description, unlike for the superfluid, the corresponding mode 
is also a ghost, in the sense that its time derivative term has the wrong sign in the action, 
see eqs. (82, 100). This makes sense since duality corresponds to a canonical transformation 
that exchanges p —» —q and q —>• p for some of the variables 15 , so that if in one description 
the potential is negative in the dual description it is the kinetic term to be negative. As 
a final remark, the quantum cut-off A of the fluid theory goes to zero like p 1 / 4 rs-/ B 1 / 8 as 
B —> 0. Therefore the fluid is formally infinitely strongly coupled at the ghost-condensate 
point 16 while the superfluid has a finite UV cut-off, and thus represents a well to do effective 
field theory. This, and the presence of two distinct phases with opposite energy density, 
represent the novelty of the ghost condensate model with respect to ordinary fluids. 

13 By eq. (102) we have p = B 3 F f' + = G(X) — G(c*) > 0. 

14 This follows from Euler equations where acceleration along fluid moton is proportional to Vp/p. 

15 Example B and E in electric-magnetic duality. 

16 For ordinary fluids associated to weakly coupled theories, the physical cut-off A p h ys where the hydrody¬ 
namic description breaks down is always smaller than p 1 / 4 and is determined by the inverse mean free path. 
For instance in a hot plasma we have A p h ys ~ g 2 T <C T ~ p 1 / 4 . 
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6 Degenerate cosmic supersolids 


Throughout the paper we have been focussing on the scalar field excitations, neglecting the 
backreaction of gravity. As we explained in the beginning, this is consistent as long as we limit 
ourselves to distances r <C rc = Mpi/ A 2 , where A is the characteristic energy scale of the 
scalar Lagrangian. The gravitational length rc normally controls both the mixing between 
the Goldstone boson and the graviton, basically the graviton mass, and the curvature of space 
time 17 . Because of the coincidence of these two length scales, the mixing between graviton 
and Goldstones cannot normally be discussed around flat space, and is thus qualitatively 
different from the Higgs phenomenon in ordinary gauge theories. One can however think of 
peculiar situations where T /tl , -C A 4 can be practically neglected, while the graviton mass 
is still of order 1 /rc- In this situation, the propagation of gravity is modified at distances 
tq where curvature can be neglected, thus defining infrared modifications of gravity that 
are the analogue of the ordinary Higgs mechanism. These theories include the Fierz-Pauli 
model, the ghost condensate model [11] and the general class of massive gravities of [13] 
and are indeed particular cases of the class of scalar theories that we have been focussing 
on. In these models, the difficulties and peculiarities associated to the dispersion relation of 
the Goldstone fields are also directly connected to the difficulty in violating NEC in general. 
Indeed when T /xu = 0 eq. (52) simply reduces to L(p,q) — L(q,p) = 0. Since the matrix A 7 
is non-zero, it immediately follows that either T_ or T + , or both, do not exist. Therefore, 
these theories can only be acceptable in limiting situations in which some dispersion relation 
is degenerate, i.e. either u = 0 or k = 0, and the corresponding mode does not propagate. 
This is in agreement with the results of Ref. [13] where the conditions for stability were 
studied in detail. Because of these degenerate dispersion relations, the stability of these 
models is threatened both by higher order effects, ex. higher derivative terms, and by slight 
deformations of the background. However in Ref. [13] it was also shown that in some cases the 
singular dispersion relations can be protected by internal symmetries of the scalar manifold, 
thus ensuring in a natural way the stability of the system. 

Now, in the presence of these non-propagating modes, the general discussion of the 
stability conditions in section 4 does not apply directly, since there we had assumed that the 
determinant of the kinetic matrix Ljj(p)/k 2 had 2N single roots. Therefore, as in [13], we 
must take a case-by-case approach, and do a very similar analysis to what was done there, to 
look at what happen if we relax the assumption that the energy-momentum tensor vanish. 
In particular we would like to know if there are stable cases when w < — 1. 

The case we study is not as general as the one considered in section 4 as we consider only 
4 fields (j)j with / = 0,1, 2, 3 in a rotationally invariant situation. However, the assumption of 
rotational symmetry is not a limitation when thinking of cosmological applications, while the 

17 We are treating the function F (see eq. (3)), as well as its derivatives and d(/)/A 2 as 0(1) quantities. 
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4-fields corresponds to the simplest case where the full translation group is broken. Setting 
as before A = 1 the background values for 0/ are taken to be 

0o = cut (107) 

0j = (3x l i = 1, 2, 3 . 

Here, a and (3 are arbitrary. By the discussion in section 5 this system consists of a superfluid 
(0o oc t) interacting with a solid (0j oc x l ), hence the denomination supersolid. 

The Lagrangian for the fluctuations 7To, 7q around the background (107) can be written, 
in the notation we have been using so far, as 

(70,o + 270) OOO cr) (d t 7io) 2 + Foo + F 0i0i ^-^j (<9j7r°) + ^ a 2 + —^ (<07q) 2 

- ( + #0i,0i) •' + ( -0 " + y 3f '*’ J:l 48 ) (ftl Ujf 

+ <i 1 ( 3f %~ faJi )w ! (108) 

We can also calculate the energy-momentum tensor in this background, 

= crj^v + 2 7000,0 0,0 + ■^F ii p 2 5 fl j5 U j^j (109) 

The null energy condition is violated if 

P + P — ot 2 70)o + g70j/3 2 < 0 (HO) 

To analyze this system, it is convenient to rescale 7To and 7Ti and rewrite the Lagrangian as 
in [12, 13], 

2 m 2 0 ( d t ir 0 ) 2 + m\ (d t i q) 2 + [m\ - (p + p)] (^7r 0 ) 2 + (4m 2 - 2m 2 ) 7T 0 <0<07T; 

- ml ( di-Kj f - (m 2 - 2m 2 ) (<07r.j) 2 (111) 

The only difference with the Lagrangian of [13] is the appearance of p + p in the coefficient 
of (Sj7To) 2 . As we explained before, for rotational invariant systems the positivity of the 
Hamiltonian is a necessary and sufficient condition for stability. In eq. (Ill) for p + p < 0, 
the gradient energy is made positive by having m 2 sufficiently large and negative. In this 
case, however, the coefficient of (<07r*) 2 would have the wrong sign, and the Hamiltonian 
would still be not positive definite. We thus recover our general result that, in the rotational 
invariant case, we cannot achieve stability with all the modes propagating. We have thus to 
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consider the degenerate cases. The analysis is immediate in the vector sector (the modes nj 
that satisfy d l nj = 0) where the Lagrangian is simply 


-^vector = ml (do7rfy - ml (dinj)' 


In the scalar sector, consisting of the modes n 0 and n L defined by n, 
Lagrangian can be written in momentum space as 


k 2 (vr 0 


^l) 


/'mlv 2 + rn\ — (p + p) 
\ i l-i is 


—XfJLlS 
2 2 

m\v — r] 



( 112 ) 

l/y/ff<9,:7r L , the 


(113) 


with is — (jj/k, /i — (2 m\ — m\) and r / = 2 (ml — m|). The analysis of the stability of 
this theory is very similar to the analysis done in [13]. The absence of ghosts and gradient 
instabilities constrains the sign of the coefficients in the diagonal entries of the kinetic matrix 


mg >0 ml > 0 , —ml + ijp + p) > 0 , rj > 0 


(114) 


This condition agrees with our general result: for p + p < 0 there cannot be stability if all 
modes are normal. We must however analyze with care the situation where some of the 
modes become degenerate, corresponding to some root v going to either 0 or oo. If the 
unstable modes can be forced to these limiting dispersion relation by consistent symmetry 
requirements, then the system has a chance of becoming acceptable. The limiting situations 
arise when one or more of the coefficients in eq. (114) equal 0. For instance for ml = 0 one 
of the roots is v = oo and the corresponding held has a pure gradient quadratic action. By 
the arguments in refs. [12, 13], we may thus forget about this mode, provided it does not 
reappear when higher derivative terms are considered or the solution is slightly perturbed. 
The other mode has a normal dispersion relation and is not a ghost as long as m\, rj > 0, 
irrespective of p + p < 0. However we know of no symmetry that can naturally enforce 
777g = 0 without introducing extra unprotected degenerate modes, so this possibility does 
not seem viable. Similarly, the case rj = 2{ml — ml) = 0 allows to send the unwanted root at 
u = 0 while keeping p+p < 0. But also this second choice cannot be enforced by symmetries. 

In Ref. [13] several internal symmetries of the scalar held manifold which protect the 
required degeneracies at T /1U = 0 were presented. These symmetries can be thought of 
as generalizations of the volume preserving diffeomorphism invariance which is relevant for 
ordinary fluids. In what follows we would like to survey the various symmetries to see 
whether they can make the p + p < 0 case acceptable. 

First we consider the symmetry 

0^^+e(0°), (H5) 


which forces mi = 0. This symmetry implies the presence of three instantaneous sound waves 
with dispersion relation k = 0, so we may call this system a rigid supersolid. Notice that 
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while for a fluid the internal symmetry (see eq. (71)), is making the solid totally deformable, 
eq. (115) is making the solid infinitely rigid. Cosmological and phenomenological properties 
of rigid supersolids with zero energy-momentum were studied in [16, 28]. In general symmetry 
(115) implies that the vector modes have a dispersion relation of the form k = 0. The 
Lagrangian for the scalar modes can be written as 


k 2 (vro 



f 2 m 2 0 - (p + p) 
\ 2 im\v 


-2im\v \ /7T 0 \ 

2 (ml - mi)J [ttJ ’ 


(116) 


According to the previous discussion one mode has dispersion relation k 
one 


UJ 


2 


P + P 



k 2 . 


0 and the second 
(117) 


For real u the sign of the corresponding residue is the same of p + p so that this mode is a 
ghost when NEC is violated. 

Another interesting symmetry that works in case of massive gravity is 


f - <W) 


(118) 


that is the full group of diffeomorphisms of the 0* submanifold. Notice that, although this 
is a very big symmetry, one can still write a non-trivial action thanks to the presence of 0 °. 
This symmetry implies m 2 = m 3 = m 4 = 0, and the presence of three sound waves with 
zero propagation velocity, u> = 0. Consequently, a system possessing such a symmetry may 
be called a superdust. Both vector modes in the superdust have dispersion relation u = 0. 
The scalar sector is characterized by the following kinetic term, 


k 2 (mo 


*1) 


(2mlu 2 + m 2 — (p + p) 
\ — im\i / 


im\v\ 

( % 

m\u 2 ) 

w 


(119) 


The eigenvalues of the kinetic matrix have zeros at to = 0 and, uj 2 = k 2 (p + p)/2ttTq. Again, 
in the case of p + p < 0, gradient instabilities can be avoided by taking mi < 0, but then the 
propagating mode is a ghost. Similarly, the symmetries 0° — > 0°+£ o (0°) and 0° — > 0°+£°(0*) 
do not lead to stable situations. 

We conclude that for p + p < 0 not even the extra symmetries that helped make sense of 
modifications of gravity can help avoid instability, p + p = 0 emerges once again as a truly 
limiting case even for this class of generalized cosmic solids. 


7 Towards w < — 1 

One of the motivations to study NEC violations in effective field theories is to understand 
whether the deceleration parameter of the cosmological expansion can be strongly negative, 
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w < — 1. It is then natural to ask what are in this respect the implications of our results. Is 
there still room for a cosmology with increasing 18 Hubble parameter? 

Our results make it very unlikely that such a cosmology may be driven by some NEC 
violating held theory. However one should keep in mind that our considerations, strictly 
speaking, apply only to scalar theories in the leading derivative approximation. For instance 
we did not consider situations where higher derivative terms in the effective action are 
important. An example of a non-trivial higher derivative theory from which to start this 
study could be offered by the dynamics of the brane bending mode in the DGP model [6, 
15]. Moreover our analysis of the symmetries which could protect the degenerate dispersion 
relations did not cover the most general case, even for scalar theories. Also it may be 
possible to construct models where instabilities are present, but slow enough to be benign. 
On the other hand, we expect our results to be quite robust, especially because they relate 
NEC violation to the instabilitites in a sector of the theory, the Goldstones of spontaneously 
broken space-time translations, that exists whenever there is a non-trivial held background. 

Here we would like to discuss yet another possibility for having w < — 1. The basic 
idea is that the accelerated expansion of the Universe is due to a genuine “modification of 
gravity”, by which the observed metric depends on new degrees of freedom and thus does 
not satisfy the Friedmann equation. For instance, one may imagine a probe 3-brane moving 
in a curved bulk space. This brane may be accelerated by forces of non-gravitational origin 
(e.g., it may be charged under some form fields). Then the induced metric on the brane will 
evolve in time even in the absence of an energy-momentum tensor localized on the brane. 
Clearly, the challenge for proposals of this type is to explain why gravity is described by the 
4-dimensional Einstein theory at short length scales. Scalar tensor theories of gravity are 
also simple enough arena to toy with this idea. This was indeed done in Ref. [30] (for earlier 
work see [31]) with the conclusion that “only highly contrived models would lead observers 
to measure w < — 1”. So our first purpose here is to understand in simple terms the basic 
reason for contrivedness in such models. This will allow us to propose a relatively simple 
candidate model where w < — 1 can be achieved without fine-tuning 19 . A minimal 4d setup 
possessing all the ingredients needed to implement this idea is described by the following 
action 

S[g, ip, 0] = Mpj J d 4 xy/=gR + S m [ip, f((p/M)g llv ] + A 0 [0, g^] . (120) 

18 Notice that there are proposals [8, 29] that mimic the presence of a dark energy component with w < — 1 
while the Hubble parameter still always decreases with time. In other words in these models the total (p + p), 
dark energy plus dark matter, is still positive. 

19 More precisely, in the model we shall present there is no extra fine-tuning associated to w < — 1. Still 
there is a necessity to fine-tune the potential of certain scalar held to make it very smooth. This is a 
consequence of the presence of a dynamical “dark energy” with w ^ — 1. This type of fine-tuning is inherent 
in most “quintessence” scenarios. 
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Here 0 collectively denotes the Standard Model fields, while 0 is an extra scalar field. The 
goal is to construct an action of the form (120) allowing for a solution with the following 
properties. 

(A) The field 0 changes in time and the physical metric /(0)<7 M „ has the FRW form. 
The main source of visible FRW expansion is the time dependence of the factor /(0/M). In 
particular, this implies that the energy density of the field 0 itself should be small enough 
not to cause significant backreaction. 

(B) Deviations from Einstein theory at short length scales due to coupling of the field 0 
to the matter energy-momentum tensor are small enough. 

In other words, the field 0 is a non-gravitating cosmic clock. Let us see that a conventional 
scalar field with the action 

S<P = j d 4 x^g(d^(j)) 2 - V(<j>) 


is very unlikely to play such a role. In order to proceed, let us first assume that the interaction 
of 0 with matter is weaker than gravity, M 3> Mpp In this case, condition (B) is satisfied 
and the time evolution of the /(0/M) factor in the physical metric gives a contribution to 
the observed expansion rate of the Universe 



(m) 


where we assumed that f(x) is a generic smooth function of order one. On the other hand 
there is a contribution to the expansion rate of the Universe from the energy-momentum 
tensor of the field 0 itself, which is at least of order 



( 122 ) 


The latter contribution gives w > — 1 and dominates the first one in the considered regime, 
M Ad pi. 

Let us now discuss the case in which the clock field 0 is coupled to matter more strongly 
than gravity, M < Mpp In this case measurements of the deflection of light by the Sun 
combined with tests of the Newton law at sub-millimeter scales imply that 0 should have 
rather big mass, 

77 i0 >0.1 mm -1 . (123) 

Now, during a Hubble time 0 varies by 


A0 ~ M 


(124) 
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so that the potential V(</>) gives a contribution to the energy density of order 

p'P ~ mjA</> 2 ~ m\M 2 . (125) 

In order for this energy to be subdominant in today’s Universe one needs 

M < ■ (0.1 mm) -2 . (126) 

In other words, the coupling of the clock held to matter becomes non-perturbative below 
0.1 mm. It is hard to imagine how such a theory could make sense. 

To summarize, we see that in both cases the essence of the problem is that the gravi¬ 
tational backreaction of the cosmic clock held (j) tends to give a larger contribution to the 
expansion rate than the direct coupling 

One way to avoid this problem might perhaps be to consider the case M < Mpi and try 
to hide the held <fi through the chameleon mechanism [32], instead of introducing the large 
mass m 

Another possibility, instead, is to consider a clock held coupled to matter more weakly 
than gravity, M 3> Mpi, and to try to minimize its gravitational backreaction. A natural 
candidate to the role of non-gravitating clock is the ghost condensate 

SM9^} = f d 4 x^A 4 G(X) (127) 

where X = (<9 M 0) 2 /A 4 as in sect. 5.3 but we also assume 

G\X 0 ) = G(X 0 ) = 0 . (128) 

for some X 0 . Without loss of generality we may set X 0 = 1. Note that the direct coupling to 
matter breaks the shift symmetry of the ghost condensate. Then, in general, matter loops 
will generate a potential for the held 0, which is too large for the construction below to work. 
The hne-tuning needed to cancel this potential is not related to the requirement w < — 1 
and is similar to one present in most quintessence scenarios. Then one finds the following 
“cosmological” solution, 

4 > = X 2 t ( 129 ) 

with matter helds in their vacuum state. This solution provides an extreme implementation 
of the original idea, in which the dependence on the held (j) is the only source of time evolution 
for the physical metric. The Hubble parameter is given by (121) and is of order 

H f ~ A 2 /M . (130) 

Its time evolution can be fairly arbitrary, in particular it can grow with time. To check that 
this model may present interest for cosmological applications let us check that values of H of 
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the order of the current Hubble rate Hq can be achieved without contradicting to the tests 
of general relativity. Two type of effects should be suppressed in this model. First there is 
a modification of gravity dne to the mixing between the metric and the ghost condensate 
mode. Second there is an extra scalar force arising from the direct coupling of the trace of 
the energy-momentum tensor to the ghost condensate held. Mixing with gravity gives rise 
to a Jeans-like instability with characteristic time-scale 

to ~ MpJA 3 . 


There are no sizable effects dne to this mixing as long as t c is larger than the current age 
of the Universe (and even for somewhat shorter values of t c ). This implies that as large 
values as A ~ 100 MeV are allowed. Then eq. (130) implies that the parameter M which 
determines the strength of the direct coupling of the ghost condensate held to matter can be 
well above Mpi still providing a realistic value of H. The only subtlety here is that in the 
limit in which gravity is decoupled the ghost condensate has the unusual dispersion relation 

a, 2 - U = 0 , (131) 


with the /c 4 -term coming from higher-derivative terms. Naively, this dispersion relation leads 
to the puzzling conclusion that the exchange of the ghost condensate held leads to a static 
potential which grows proportionally to the distance from the source. This would imply the 
breakdown of perturbation theory far from the source even for very large values of M. A 
resolution to this problem is that the dispersion relation (131) implies that the velocity of 
the scalar quanta 


(132) 


is very low for small values of k. As a result a very long time ~ L 2 A is needed for this 
potential to be established at a large distance scale L. More precisely accounting for the 
mixing with gravity the dispersion relation (131) is modified to 


k A k 2 A 2 


w ■ F + Mi 


= 0 


PI 


(133) 


The last term gives rise to the Jeans-like instability discussed above and dominates over the 
second one for 


k < m = 


A 2 


M, 


pi 


(134) 


Consequently, for time-scales shorter than t c the linear potential is established only for 
distances shorter than m _1 . At larger distances mixing with gravity trades this growing-with- 
distance-potential with an instability that develops in time. For small enough A < 100 MeV 


and large enough M 3> Mpi both effects can be neglected. 
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Note, that in the presence of massive matter sources non-linear effects dominate the 
dynamics of the ghost condensate [27]. These effects may introduce local spatial gradients 
of the held 0 of order di(f) ~ A 2 . Then direct coupling of the ghost condensate held to matter 
may result in local anomalous accelerations of order A 2 /M ~ H 0 . This may be an interesting 
or deadly feature of this model, depending on whether such effects are present inside the 
Solar system. The detailed study of these effects is beyond the scope of this paper. 

8 Summary 

We have studied the connection between the null energy condition and the microscopic 
stability of the system that sources the energy momentum tensor. We have not considered 
extra pathologies that may arise through the gravitational backreaction. Our attention has 
been limited to a general system of derivatively coupled scalar helds, but we expect that our 
arguments can be generalized to include vector helds. However, we have explained that the 
class of theories we consider includes relativistic solids and fluids, and is thus relevant to 
cosmology. 

We have hrst worked under the assumptions that all gapless scalar modes have a non¬ 
degenerate dispersion relation ui = vk with v ^ 0, oo. Our main result then is that a 
violation of the null energy condition does not strictly imply instabilities (ghosts or imagi¬ 
nary frequencies) in the system. We have proven that by producing counter-examples, that 
is acceptable effective held theories with a positive definite Hamiltonian for quadratic pertur¬ 
bations around a background whose energy momentum violates the null energy condition. 
A necessary feature of all the counter-examples is the anisotropy of the background and, 
perhaps more importantly, the presence of superluminal modes. In fact we proved that for 
systems that are either isotropic or do not feature superluminality, a violation of the null 
energy condition always implies an unescapable instability. From a purely effective held 
theoretic viewpoint, in the presence of a Lorentz violating background, there does not seem 
to exist any limitation on the presence of superluminal modes. Indeed many (most?) of 
the models that have been considered in recent years in connection to cosmic acceleration 
[33, 6, 11] feature superluminal excitations around some background. Our study suggests 
that effective held theories that allow superluminality are subtly distinct from those which 
do not. The apparent paradoxes that we discussed in connection to causality are another 
sign of this distinction and a simple, but deep, characterization of it will be given in Ref. [34], 

We have then considered isotropic systems in which some mode may have a degenerate 
dispersion relation, either u = 0 or k = 0. We expect such systems to be acceptable only 
when the degenerate dispersion relations are protected by symmetries. We have studied 
various internal symmetries of the scalar held manifold. These symmetries nicely generalize 
to relativistic solids the internal symmetry that distinguishes ordinary liquids from ordinary 
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solids (see paragraph below). Although we have not been able to perform a completely 
general analysis, we believe we have given fairly convincing evidence that a violation of the 
null energy condition leads to instability even in the presence of such symmetries. It is 
interesting, in this respect, that the so called theories of massive gravity, broadly defined by 
having a non-trivial field background while the energy momentum tensor vanishes exactly, 
lie right at the boundary between theories that violate and theories that respect the null 
energy condition. This feature largely explains the difficulty to make such theories viable. 

Aside the above results, two other directions were partially investigated. In sect.5 we 
considered the subclass of scalar field theories involving three scalar fields (j) 1 , / = 1,2,3, 
with space coordinate dependent expectation values (0 J ( x )) = x 1 . We explained that these 
systems are naturally interpreted as solids filling the universe. We characterized the limit 
where the solid becomes a fluid by the appearance of a volume preserving reparametrization 
invariance of the scalar field manifold. We then generalized to the relativistic case in a 
Lagrangian formulation some of the well known results of fluid mechanics, like Kelvin’s 
theorem. We also discussed a sort of duality relation between the recently proposed ghost 
condensate model and an ordinary diluted fluid. In this respect the novelty of the ghost 
condensate models partially lies in its being a weakly coupled theory, which is not the case 
for the corresponding diluted fluid model. In our opinion, the results in sect. 5, while not 
completely original and not fully clear cut, are still useful as they put in a wider perspective 
the classes of scalar theories we have been working on. 

Finally in section 7 we briefly stepped into another class of theories, where one scalar with 
non-trivial time dependent background is directly coupled to the matter energy momentum 
and effectively mimics a cosmic acceleration with H > 0. We have simply explained the two 
difficulties that have made past attempts of this type very clumsy. These are due to the tests 
of general relativity and to the gravitational backreaction of the metric on the scalar field 
energy momentum. The two difficulties point to a simple solution, where the scalar field is 
a ghost condensate field, for which the energy momentum can be exactly zero thus causing 
no backreaction. The possible difficulties of such a model lie in the non-linear dynamics 
of the ghost condensate in the presence of massive bodies, though it has reasonably been 
argued that such non-linearities may not necessarily lead to trouble. Our simple proposal 
potentially opens the way to genuine modifications of gravity that make w < — 1 consistent, 
but further study is needed. 
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A Stability of the NEC-violating systems of section 4.1 

In section 4.1, we presented two Lagrangians which were claimed to be stable despite violating 
the NEC. Here we sketch the stability analysis of these systems. 

A.l Two dimensional example 

The Lagrangian for the two dimensional example is given by 

C = (h 00 + B 11 + 2(1 + a 2 ) B m - (a 2 + b 2 ) B 00 B m + (a 2 + c 2 ) . (135) 

with 

b 2 < 2, c 2 < 2 . 

To study the stability of the fluctuations, it is useful to know the expansion of the meson 
field B IJ , 

B 00 = 2d 0 n° + d^d^n 0 
B 11 = —2 dm 1 + d^d^n 1 
B 01 = (don 1 - dm 0 ) + dyd^ 1 

The Lagrangian for the fluctuations can then be written as 

(A7r 0 ) 2 + (^tt 1 ) 2 + 2(1 + a 2 ) («9 m 7t°) (d^n 1 ) 

— (a 2 + b 2 )2d 0 n 0 d 0 7r 1 + (a 2 + (?)dm l dm° + (a 2 + c 2 )2<9 0 7r°(9i7r 0 — (a 2 + c 2 )dm 1 do'K 1 . 

(136) 

The first line of (136) comes from the first three terms of (135) which are linear in B 1J . Those 
are the only term of (135) that contribute to the background energy-momentum tensor and 
are responsible for the NEC violation. We can see that taken alone, they would form an 
unstable system as both the kinetic and gradient matrices of the first line of (136) have 
negative eigenvalues, due to the mixing between 7r° and ti 1 . But this defect is cured by the 
terms B 00 B 01 and B ll B 01 that provide terms in the fluctuation Lagrangian that suppress 
the kinetic and gradient mixing respectively. They also give terms with one space and one 
time derivative which do not affect the positivity of the Hamiltonian (cf. eq. (25)). 

A.2 Three dimensional example 

The Lagrangian written in terms of B 1J has form 

C = T u B IJ + (a 2 [(H 01 ) 2 - (H 11 ) 2 - (H 12 ) 2 ] + ch sh (1 + a 2 ) [B 02 (-B 11 + B 22 ) + 2 H 01 H 12 ]) 

(137) 
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with a 2 ~ b 2 ~ a 2 ~ 1 and Tjj as given by eq. (64). As in the two dimensional example, 
terms linear in B IJ give a Lagrangian for the fluctuation of the form Tjjd^m 1 d^rr 7 that is 
Lorentz invariant and unstable, due to the — 6 2 (tt 1 ) 2 that it contains. And as before, this 
can be cured by adding the terms that are quadratic in B IJ . To understand what these 
terms do, it is useful to keep in mind the expansions of the B 1J to linear order in it 1 , 


B 00 = 2 d 0 n° 

B n = -2 dm 1 
B 22 = -2 dm 2 
B 01 = dm 1 — dm° 

B 12 = —dm 2 — dm 1 
B 02 = d 0 TT 2 — dm 0 


The kinetic part of the Lagrangian for the fluctuations is given by 

K = (ch 2 a 2 + sh 2 ) (d 0 n 0 ) 2 +2chsh(l+a 2 )d 0 n 0 d 0 n 2 +(ch 2 +a 2 sh 2 ) (<9 0 7r 2 ) 2 + (a; 2 — b 2 ) (dovr 1 )^ , 

(138) 

while the gradient matrix is 


G = (ch 2 a 2 + sh 2 ) (dm 0 ) 2 + ( sh 2 a 2 + ch 2 ) (dm 2 ) 2 + (a 2 — b 2 ) (dm 1 ) 2 

— a 2 ((9i7t 0 ) 2 + dm 1 dm 2 — (dm 2 ) 2 (139) 

The role of the term a 2 (B 01 ) 2 in (137) is to make the kinetic matrix positive definite, 
which is the case if 

a 2 -b 2 > 0 . (140) 

The term ch sh (1+a 2 ) [B 02 (—B u +B 22 )+2B 01 B 12 ] removes mixing between it 2 and 7r° in the 
gradient matrix. This is needed because with the mixing the gradient matrix, despite having 
entries of order chsh , has an eigenvalue that is small and the term — a 2 (<9i7r 0 ) 2 contained in 
a 2 (B 01 ) 2 could make this matrix non-positive definite. Note, that mixing between n 2 and 7r° 
is still present in the kinetic matrix. As a result this matrix has a small eigenvalue leading 
to the highly superluminal mode in the spectrum in agreement with our results. 

Finally, the terms —a 2 ( B n ) 2 — (B 12 ) 2 are there to change the sign of the gradient term 
for (dm 1 ) 2 ■ The terms on the last line of (139) are of order 1 and subdominant compared 
to the terms of order ch, sh and do not affect the positivity of the gradient matrix. 
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